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TRANSVERSE MOMENTUM BROADENING DUE TO THE 
MULTIPLE SCATTERING 

XIAOFENG GUO 

Depariment of Physics, Columbia University, New York, NY 10027, USA 

Using the Drell-Yan process in hadron-nucleus collisions and deeply inelastic 
lepton-nucleus scattering (DIS) as examples, I show that the transverse momentum 
broadening can be expressed in terms of four-parton correlation functions. I argue 
that jet broadening in DIS provide an excellent measurement of the four-parton 
correlation functions and a test of QCD treatment of the multiple scattering. 


1 Introduction 

Most high energy collisions can be described by a single hard scattering. By 
studying the single scattering, we can extract the parton distribution func¬ 
tions and test the short-distance dynamics of strong interactions. The normal 
parton distributions have the interpretation of the probability distributions to 
find a parton within a hadron. On the other hand, in order to study quan¬ 
tum correlations of multi-partons inside a hadron, we need to investigate the 
multiple scattering. The multiple scattering is a high twist effect and is small 
compare to the single scattering. However, in collisions involve nucleus, the 
nucleus environment will enhance the multiple scattering effect. For example, 
when a parton propagates through a nuclear matter, the average transverse 
momentum is broadened due to the multiple scattering. Nuclear dependence 
of such broadening provides an excellent probe to study QCD dynamics beyond 
the single hard-scattering picture. 

A reliable calculation of the multiple scattering in QCD perturbation the¬ 
ory requires to extend the factorization theoremlH beyond the leading power. 
Qiu and Sterman showed that the factorization theorem for hadrorthadron 
scattering holds at the first-nonleading power in momentum transfera, which 
is enough to study double scattering processes in hadronic collisions. Ap¬ 
plying this generalized factorization theorem, Luo, Qiu and Sterman (LQS) 
developed a consistent treatment of the multiple scattering at the partonic 
levela. LQS expressed the nuclear dependence of di-jet momentum imbalance 
in photo-nucleus collisions in terms of four-parton (twist-4) nuclear correlation 
functionsQ. In the following, I will use the the deeply inelastic lepton-nucleus 
scattering (DIS) and the Drell-Yan process in hadron-nucleus collisions as ex¬ 
amples to show that transverse momentupi broadening is proportional to the 
similar four-parton correlation functions 0. 

The predictions of the multiple scattering effects rely on accurate infor- 
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mation of these four-parton correlation functions. The four-parton (twist-4) 
correlation functions are as fundamental as the normal twist -2 parton distribu¬ 
tions. Just like normal parton distributions, multi-parton correlation functions 
are non-perturbative and universal. QCD perturbation theory cannot provide 
the absolute prediction of these correlation functions. They can be measured in 
some processes and be tested in other processes. Using the Fermilab E683 data, 
LQS estimated the size of the relevant twist-4 parton correlation functions to 
be of the order of 0.05 — 0.1 GeV^ typical twist-2 parton distributions 

□. However, Fermilab and CERN dataQ’0 on nuclear enhancement of the aver¬ 
age Drell-Yan transverse momentum, A(g|,), prefer a much smaller size of the 
four-parton correlation functions □’□. This discrepancy may result from differ¬ 
ent higher order contribution to nuclear enhancement of the average Drell-Yan 
transverse momentum Q’tu and the di-jet momentum imbalance. It is neces¬ 
sary to study the higher order corrections to these two observables in order to 
test QCD treatment of multiple scattering. Meanwhile, it is also important to 
use jet broadening in_DIS as an independent measurement of the four-parton 
correlation functionsu. 

In the following, I first present the derivation of the transverse momen¬ 
tum broadening in Drell-Yan process and jet broadening in DIS. I then argue 
that jet broadening in DIS is an excellent observable to study the four-parton 
correlation function. 


2 Transverse momentum broadening for the Drell-Yan pairs 

Consider the Drell-Yan process in hadron-nucleus collisions, h{p') + A{p) —» 
(g) -|- X, where q is the four-momentum for the virtual photon 7 * which 
decays into the lepton pair, p' is the momentum for the incoming hadron and 
p is the momentum per nucleon for the nucleus with the atomic number A. 
In order to extract the effect due to the multiple scattering, we define the 
Drell-Yan transverse momentum broadening as 

A(g2 ) ^ - {qlt^ , (1) 


with qx the transverse momentum of the Drell-Yan pair, and the the averaged 
transverse momentum square is defined as 


(Qt) 


2 \hA 



dOhA 

dQ^dq^ 


dcThA 

dQ^ 


( 2 ) 


In Eq. (Hj), Q is the total invariant mass of the lepton pair with = q^. From 
our definition, A(g|.) represents a measurement of QCD dynamics beyond the 
traditional single-hard scattering picture. The nuclear dependence of A{q^), 
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defined in Eq. (|^) , is a result of the multiple scattering between the parton from 
incoming beam and the nuclear matter before the Drell-Yan pair are produced. 

If we keep only the double scattering contribution and neglect contribution 
from the higher multiple scattering, we have 


MQt) 



d^hA 

dQ'^dq^ 


/ 


dan A 


(3) 


where superscript “D” indicates the double scattering contribution. Fig. || 
shows the leading order double scattering between the parton “f” from incom¬ 
ing beam and the nucleus. According to the generalized factorization theorem 
□, the double scattering cross section can be expressed as 


d^hA^t+i- 



(4) 


with ddfj^ the parton level double scattering cross section, and 
1 


= 2x', 


dxdxidx^ j dld^T^d^x,Xi,X2, kx) 


xH{x,Xi,X2,kT,x'p',p,q) . 


(5) 


In Eq. (Q), cj)f/fi{x') is the parton distribution from the hadron. In Eq. (||), H 
is the hard partonic part, and T is the hadronic matrix element: 


T^^\x,Xi,X2,kT) 

— f dpi dy 2 d Ut ixp+y- ixip+(v 7 -y~) ix2.p+v~ ikT-VT 

J 4tt 27r 27r (2^)2® 

X {pA\A+{y^,OT)ipq{0)j^i’q{y~)A+{y:[,yiT)\PA) ■ (6) 

Because of the exponential factors in Eq. i), the position space integration, 
dy~’s, cannot give large dependence on the nuclear size unless the parton 
momentum fraction in one of the exponentials vanishes. If the exponential 
vanishes, the corresponding position space integration can be extended to the 
size of the whole nucleus. Therefore, in order to get large nuclear enhancement 
or jet broadening, we need to consider only Feynman diagrams that can provi^ 
poles which set parton momentum fractions on the exponentials to be zeroBu. 
Other diagrams that do not provide such poles will be suppressed by a large off- 
shell propagator and hence are not leading contributions. At the leading order 
in as, only diagrams shown in Fig. ^ have the necessary poles. These diagrams 
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Figure 1: Lowest order double scattering contribution to the nuclear enhancement of Drell- 
Yan (g^); (a)symmetric diagram; (b) and (c): interference diagrams. 


contribute to the double scattering partonic part H(x,Xi,X 2 ,kT,x'p',p,q) in 

Eq. (D. 

For the leading order diagrams shown in Fig. |^, the corresponding partonic 
parts have two possible poles from the two propagators. For the diagram shown 
in Fig. 1^, the partonic part has the following general structure 


1.2 

xS{x + xi - ^ - ^)d{qT - k^) . (7) 

x's x's 

The two 5-functions are from the phase space. One of the 5-functions can be 
used to fix dx integration in Eq. (||), and the two poles in Eq. (|^) can be used 
to perform the contour integration for dxi dx 2 - We have 

dda (x{ATT^)e{y--y^)0{-y2)S{q^-kl)e^^'^/^'''P^y , ( 8 ) 

with T = js. The two 0-functions are the results of the contour integration. 

For the diagrams shown in Fig. |pD and|^c, we have slightly different phase 
space factors and different poles. Similarly, after integrating over parton mo- 
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mentum fractions, dxdxi dx 2 , we have 

daj^ (x{-4TT^)e{y- -y-)0{y- , (9) 

and 

da° oc (-dTT^) 6»(yf - y^) 0{-yi ) 5{q^)e' . (10) 

The spinor trace gives the same numerator for all three diagrams. There¬ 
fore, the total contribution to transverse momentum broadening is proportional 
to -I- ddf* -|- dd^ , and this sum has the following feature 


dQ"^ dq^ 


oc dd^ + ddj^ + dd^ 


« 0iy~ - 2/r)6'(-?/^) [Siq^ - fcr) - 

+ [s{y~ - yi) Q{-y2) - - yi) ^iy~ 

-^{yT -y2 )^{-yi)] Hyr) ■ 


02 ) 


( 11 ) 


It is clear from Eq. 0 ) that for the inclusive Drell-Yan cross section, da/dQ^, 
the double scattering contribution, daf;’j^/dQ^, does not have a large depen¬ 
dence on the nuclear size. The integration over dq^ eliminates the first term 
in Eq. dO) , while the second term is localized in space if t/x' is not too small. 
When the t/x' is finite, and is large, exp[i (r/a;')p+j/“] effectively restricts 
y~ ~ 1/((t/x')p+) —> 0. When y~ 0, the combination of the three pairs of 
d-functions in Eq. & vanishes, 


dc^hA 

dQ^ 


= J dq1^(^ 


dc^hA 

dQ'^dq^ 


oc [diy - 01 ) 6'(-j/2 ) - ^{,02 - 01 ) 0 {y~ 
^0 as > 0 . 


-y2)-d{yi -y2)d{-yi)\ 

( 12 ) 


Physically, Eq. ( 0 ) says that all integrations of y ’s are localized. Actually, 
at the leading order, the term proportional to 9{y~ — yf)d(—y^) — 6{yl// — 
y/) 0{y~ — y^) ~ d{yi — y^) d{—yi) is the eikonal contribution to make the 
normal twist-2 quark distribution for single scattering gauge invariant. Eq. (12) 
is a good example to demonstrate that the double scattering does not give a 
large nuclear size effect to the total inclusive cross section. 

On the other hand, from Eq. (0), the double scattering contribution to 
the averaged transverse momentum square can give a large nuclear size effect 



f da^j, \ 

\dQ^dq^ ) 
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(13) 


cx / dq^ Qt [^{Qt ~ ^t) ~ ^{Qt)] 


Actually, k'^ in Eq. needs to be integrated first. But Eq. (0) already 
demonstrates that ^{q^) is proportional k^, which is the kick of the transverse 
momentum the parton received from the additional scattering. The bigger the 
nuclear size, the bigger the effective k^. As shown below, ^{q^) is proportional 
to the nuclear size. 

After working out the algebra, we obtain the Drell-Yan transverse momen¬ 
tum broadening at the leading order in as 


A{ql) = 


An^as \ Eg dx' ^gjhix') T^j\{Tlx')lx' 

Eg e'^q I dx' (j)glh{x')(l)g/AiTlx)lx’ 


(14) 


In Eq. (0, ipq/A^x) is the usual quark distribution inside a nucleus. T^^^j^(x) 
is the four-parton correlation function, and is given by 




27T 


dVi dy^ 
27r 


d{y -yi)d{-y2) 


X 2 iPA\FEi.y 2 )V’<?( 0 ) 7 ’^’(’<?(y )-F^“(2/i )\pa) ■ 


(15) 


The superscript (I) represents the initial state interaction, in order to distin¬ 
guish from the similar four-parton correlation function defined in Eq. j^). In 
deriving Eq. we expend S{q^ — k^) at kr = 0, known as the collinear 

expansion, and keep only the first non-vanishing term which corresponds to 
the second order derivative term. We use the factor k^k^ to convert the 
into field strength F°‘~^FP^ by partial integration. Here, we work 
in Feynman gauge. The terms associated with other components of are 
suppressed by l/p”*" compared to those with A+, because of the requirement 
of the Lorentz boost invariance for the matrix elements 0 . 

By comparing the operator definitions of these four-parton correlation 
functions and the definitions p^the normal twist-2 parton distributions, LQS 
proposed the following model! 


Tf/Aix) = , 


(16) 


where (j^f/Nix) with / = q^q,g are the normal twist-2 parton distribution 
of a nucleon, and A is a free parameter to be fixed by experimental data. 
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Applying this model, we obtain a very simple result for the Drell-Yan transverse 
momentum broadening: 


A(4> = . (17) 

Using Eq. ( 0 ) and data from E772 and NAlO on nuclear enhancement of the 
average Drell-Yan transverse momentum, A{q'^), we estimate the value of 
as 

« O.OlGeV^ , (18) 

which is at least a factor of five smaller than 0.05—0.1 GeA^, previously 

estimated by LQS from momentum imbalance of the di-jet datau. Therefore, 
it is very important to use other observables, such as jet broadening in DIS, 
to further test the four-parton correlation functions. 


3 Jet broadening in Deeply Inelastic Scattering 

In this section, we derive the leading contribution to jet broadening in DIS by 
using the same technique used to derive the Drell-Yan transverse momentum 
broadening, A(g|,), in last section. Consider the jet production in the deeply 
inelastic lepton-nucleus scattering, e(fci) -|- A{p) —> 6 (^ 2 ) -|- jet{l) + X. ki and 
^2 are the four momenta of the incoming and the outgoing leptons respectively, 
and p is the momentum per nucleon for the nucleus with the atomic number 
A. With I be the four-momentum for the jet, the averaged jet transverse 
momentum square is defined as 




/ 


dl^-l 


2 

T ■ 


d(7 eA 

dxsdQ'^dl'^ 


dd^A 

dxsdQ'^ 


(19) 


where xb = Q^/(2p-q), and q = ki — k 2 is the momentum of the virtual photon, 
and = —q^. The jet transverse momentum It depends on our choice of the 
frame. We choose the Breit frame in the following calculation. Similar to 
the Drell-Yan transverse momentum spectrum, da/dQ'^dq^, the jet transverse 
momentum spectrum, da/dxsdQ'^dl'^, is sensitive to the A^/^ type nuclear 
size effect due to the multiple scattering. On the other hand, the inclusive 
DIS cross section da/dxBdQ"^ = f dl^ da/dxBdQ'^dl^ does not have the A^/'^ 
power enhancement. Instead, it has a much weaker A-dependence, such as the 
EMC effect and the nuclear shadowing. To separate the multiple scattering 
contribution from the single scattering, we define the jet broadening as 


A{il) ^ - {ilY^ , 


( 20 ) 
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Keeping only the contribution from the double scattering, similar to Eq. 
we have 



da 


D 

eA 


dxsdQ'^dl^ 


da e A 
dxsdQ"^ 


(D> 

( 21 ) 


The leading order double scattering diagrams contributing to jet broaden¬ 
ing are given in Fig. |2[ The double scattering contribution can be expressed 
as 


da^ (X I dx dxi dx 2 


dk^T^^\x ,Xi,X2,kT)Hi 


fll^ 


{x,Xi,X2,kT,p,q,l) , (22) 


with the matrix element 


T^^'>{x,Xi,X2,kT) 

_ f dy dpi dy2 d yr ^ ixp+y~ Jxip+(y~-y~) Jx^p^y:^ „ikT-VT 

J 47r 27r 27r (27r)2 

X A+{y2,0T)A+{y:[,yiT)tpqiy~)\pA) , (23) 

where superscript (F) indicates the final-state double scattering. The matrix 
element is jeoual to the in Eq. (^ if we commute the gluon fields with 

the quark fieldsH’EJ. In Eq. (|^), is the corresponding partonic part. Here, 
we work in Feynman gauge. For the diagram shown in Figj^, the partonic 
part has the following structure 

oc S(x + Xi — xb — ——) ■ S{Ib — k^) dl^ 

2p-q 

X--- . (24) 

X — Xb + ie X + X2 — Xb — ie 

In Eq. (^J) , the (5-function is from the phase space, and the poles are from the 
propagators. 

Following the derivation of ^{q^) in last section, first, we carry out the 
integrations of the parton momentum fractions by using the 5-function and 
two poles in Eq. (p4|), 

dcr“ oc (27r)^ 6I(?/C - y~) eijj^ ) 6(1^ - k^) ■ (25) 

Similarly, we have the corresponding integrations for the interference diagram 
in Fig. ^ 


da^ (X-{2T:f B{y 2 -y^ )e{y^ -y )6 {It) 


(26) 









Figure 2: Lowest order double scattering contribution to jet broadening: (a) symmetric 
diagram; (b) and (c): interference diagrams. 


and for the diagram in Fig. ||c 

oc -(2^)2 eiy^)0{y^ - y^)5{ll) . (27) 

Combining Eqs. (|^), ( p^ ) and (|2^), we again have the same structure as that 
in Eq. &■ Therefore, like in the Drell-Yan case, we conclude that the double 
scattering does not result into any large nuclear size dependence in the inclusive 
DIS cross section. For the jet broadening defined in Eq. (13), we drop the term 
proportional to 6»(yf - y") -yi)S{yi -y~)-S{yi ), 

which is localized as the single scattering. After carrying out the algebra, we 
derive jet broadening in DIS as 


MIt) 


^ T,q elTg/Ajxs) 

V 3 y e‘^q(t>q/A{XB) ' 


( 28 ) 


where J^q sums over all quark and antiquark flavors. In Eq. (p^, (j>q/A{x) 
is the normal twist-2 quark distribution inside a nucleus, and the four-parton 
correlation function, is defined asQ 


Tq/A{x) 


^y y 

2tt 


dyi dy2 

27T 


d{yi 


y ) d{y2) 
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>^n{pA\il’q(.^)TK (.y2)F {yi)'4’q(.y )\pa) ■ (29) 


2 \r'A^\'t^q\^) 1 ^ <7 Vi/2 ) ^ \yi ) H^q\ 

With the model given in Eq. (0), we can simplify Eq. 

^47r^as 


A(4) = 




(30) 


From Eq. (|28|), we see that jet broadening in DIS is directly proportional to 
the four-parton (twist-4) correlation function, defined in Eq. (p^). It is the 

r me four-parton correlation function appeared in dijet momentum imbalance 
Tq/A{x) and which are defined in Eqs. ( p^ and (lA), respectively, 

are same if the phase spac^mtegral are symmetric, because field operators 
commute on the light-cone laO. Therefore we conclude that jet broadening 
in DIS, transverse momentum broadening for the Drell-Yan pairs, and dijet 
momentum imbalance depend on the same four-parton correlation function. 


4 Discussion and conclusions 

From the simple expression in Eq. (p^), we conclude that at the leading or¬ 
der, jet broadening in DIS has a strong scaling property, that it does not 
depend on beam energy, and xb- This scaling property of jet broaden¬ 
ing in DIS is a direct consequence of LQS model for four-parton correlation 
functions, given in Eq. (0)- However, when x n < 0.1, j/ ~ 1/{xbP^) is no 
longer localized inside an individual nucleon Ed’Ej. Therefore, terms propor¬ 
tional to e{y- - y-)6{y-) - 0{y- - y-) e{y- - y“) - 6»(?/f - y") 6»(y-), are 
no longer localized and need to be kept for the jet broadening calculation if 
Xb is small. Consequently, the Xs-scaling of jet broadening in DIS needs to 
be modihed in small xb region. In addition, Q^-dependence may be modified 
because the four-parton correlation function Tq/A{x) and the normal quark 
distribution <l>qiA^x) can have different scaling violation. In principle, all de¬ 
pendence or whole conclusion could be modified due to possible different high 
order corrections. Nevertheless, we believe that experimental measurements 
of jet broadening in DIS can provide valuable information on the strength of 
multi-parton correlations and the dynamics of the multiple scattering. 

Similarly, from Eq. ( |l7|) , we can also conclude that the Drell-Yan transverse 
momentum broadening, A{q^), has a small dmendence on beammergy and 
of the lepton pair. Data from Fermilab E772Q and CERN NAIOQ demonstrate 
weak energy dependence. It signals that the simple model by LQS for four- 
parton correlation functions is reasonable and the leading order calculation 
given here are useful. At the same time, the observed energy depe^ence 
indicates that the high order corrections to A{q^) can not be ignored Bty. 
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In addition, Eqs. (30) and ( 0 ) tell us that at the leading order, jet broad¬ 
ening in DIS and nuclear enhancement of average Drell-Yan transverse mo¬ 
mentum, A{q'^), have the same magnitude, if the averaged initial-state gluon 
interactions is equal to the corresponding final-state gluon interactions, i.e.. 


From Eq. (p^), we see that the jet broadening is directly proportional to the 
four-parton correlation function Tg/Ai^s)- Information on xs-dependence of 
the jet broadening can provide a first ever direct measurement of the functional 
form of four-parton correlation functions Tq/A{xB, ^)- In addition, by examing 
the scaling property of jet broadening, we can directly verify LQS model of four- 
parton correlation functions, given in Eq. (p^. Future experiments at HERA 
with a heavy ion beam should be able to provide much more information on 
dynamics of parton correlations. 

In summary, we have derived analytic expressions for jet broadening in DIS 
and the Drell-Yan transverse momentum broadening in terms of universal four- 
parton correlation functions. Because the dijet data (pure final-state multiple 
scattering) and the Drell-Yan data (pure initial-state multiple scattering) favor 
two different sizes of the four-parton correlation function, measurement of jet 
broadening in DIS (pure final-state multiple scattering) will provide a critical 
test of QCD dynamics of the multiple scattering. Since at the leading order, 
the Drell-Yan transverse momentum broadening, A{q^)j^Yj and jet broadening 

a -1 DIS, A(12 are independent of the four-gluon correlation function Tg^A 

, measurement of A{q^)oY and A{l^)ois provide a direct comparison between 
the initial-state multiple scattering and the final-state multiple scattering. 
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